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Abstract 

We calculate the HQET quark-gluon vertex at one loop, for arbitrary external 
momenta, in an arbitrary covariant gauge and space-time dimension. Relevant 
results and algorithms for the three-point HQET integrals are presented. We also 
show how one can obtain the HQET limit directly from QCD results for the quark- 
gluon vertex. 
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1 Introduction 



Heavy Quark Effective Theory (HQET) is an effective field theory approximating QCD 
for problems with a single heavy quark having mass m, when characteristic momenta of 
light fields are much lower than m, and there exists a 4- velocity v such that characteristic 
residual momenta k = p — mv of the heavy quark are also small. It has substantially 
improved our understanding of heavy quark physics during the last decade Methods 
of perturbative calculations in HQET are reviewed in P| . 

In this paper, we calculate the quark-gluon vertex in the leading-order HQET (l/m°) 
at one loop, for arbitrary external momenta, in an arbitrary covariant gauge, in space- 
time dimension d = 4 — 2e. This allows us to take all on-shell limits (introducing 
additional 1/e divergences) directly. The general ci-dimensional results can also be used 
for expansion around a dimension other than 4; for example, 2-dimensional HQET was 
considered in the literature in some detail. 

A one-loop calculation of the QCD quark-gluon vertex with a finite quark mass m 
has been recently completed |^ (where references to earlier partial results can be found). 
We check how the HQET result can be obtained by taking the limit m — oo in the QCD 
result. 

Let the sum of bare one-particle-irreducible vertex diagrams in HQET (Fig. 0) be 
igQt°'r'^{k,q). The "full" momenta of the incoming quark and the outgoing one are 

p = mv + k, p' = mv + k', (1) 

where v is the heavy-quark 4- velocity (f^ = 1), and k, k' are the residual momenta. The 
momentum transfer is q = p' — p = k' — k. In the HQET limit, m —>■ oo, k k' q 
(9(1). The heavy quark propagator in HQET is 

and the elementary quark-gluon vertex is igof^v^. 

In the leading-order HQET, heavy-quark propagators and vertices do not depend on 
the component of the heavy-quark momenta orthogonal to v. Therefore, T^^{k,q) does 
not depend on k^ = k — {k ■ v)v. The only vectors in the problem are v and g, and 
(see [§) 

r^(A;, q) = V,{uj, J, q^)v^ + r,(u;, a;', q^)q^, (3) 




Figure 1: HQET quark-gluon vertex 
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Figure 2: One-loop vertex diagrams 

where 

uj = k ■ V, uj' = k' ■ V (4) 

are the residual energies, and q-v = uo' —uo. The functions Ty and Tq can be reconstructed 
from the contractions 

where 

is the 3-momentum transfer squared in the v rest frame. 

At the tree level, = v^. One-loop corrections are shown in Fig. The contribution 
F^ of the diagram Fig. ^ is proportional to v^; that of Fig. ^ has both structures. 

The contraction F'^(A;, can be simplified using the identities shown in Fig. |^. Here 
a gluon line with a black triangle at the end denotes a "longitudinal gluon insertion"; 
when attached to a vertex, it means just the contraction with the incoming gluon mo- 
mentum (note that it contains no gluon propagator!). A dot near a propagator means 
that its momentum is shifted by q. The colour structures are singled out as prefactors 
in front of the propagator differences. The circular arrow in Fig. §3 shows the order of 
indices in the colour structure of the three-gluon vertex i/°'"^. Two last terms in Fig. |^ 
contain longitudinal gluon insertions again; for them, the identities of Fig. ^ can be 
recursively used. 

Application of these identities to the diagrams of Fig. ^ is shown in Fig. ^ Here the 
non-standard vertices of Fig. |^ are igot"" and g'^f'^^'^t'^v^ . This is, of course, just the one- 
loop case of the general Ward-Slavnov-Taylor identity for the quark-gluon vertex, which 



go-- 



•v/\^ •- -<-i\/\/\/\/\» 



Figure 3: Longitudinal gluon insertions into quark and gluon propagators 
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Figure 4: Ward-Slavnov-Taylor identities for the one-loop vertex 
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Figure 5: Non-standard vertices 

is discussed in in detail. The one-loop contributions to this identity are collected 
in Eqs. (2.28) and (2.29) of ||^. They can be easily associated with the diagrams shown 
in Fig. 1^. The only term which requires some explanation is the diagram involving one- 
loop ghost self-energy contribution (the first diagram in the second line of Fig. ^). Its 
connection with the last contribution on the r.h.s. of Eq. (2.29) of is shown in Fig. ^ 
This relation can be easily verified by direct calculation. Using Fig. we can avoid 
introducing the non-standard vertex shown in Fig. Qd. 

\ 
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Figure 6: Relation between ghost self-energy contributions 
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Figure 7: Heavy-quark self energy 



For the contributions of the diagrams of Fig. ^,b, we have 

Ca 



1 - 



Ca 



V 



2Cf. 

S(c<j) — S(ct;')) + (ghost terms). 



2Cf 

Here — iS(a;) is given by the one-loop self-energy diagram of Fig. |^: 

Cf gl 

2 (471)^^/2 



:2 + (rf-3)OX(^) 



(7) 



(8) 



J(^) 



TT 



d/2 



{l-V 



UJ 



i0)(/2 + i0) 



2{-2iuY-^T{3 - d)T{d/2 - 1] 



(9) 



(see 0). In what follows, we shall not explicitly write -|-iO in denominators. Here 
^ = 1 — ao, cto is the bare gauge- fixing parameter. We can see that the Yennie gauge 
(see also in 0) is of special interest, since S(ti;) is finite at ^ = —2. Moreover, if the 
generalization of Yennie gauge to an arbitrary dimension is chosen as ^ = —2/{d — 3) Q], 
then in the Abelian case (and, in particular, at one loop), the heavy-quark self energy 
vanishes (see for a tutorial). The two- loop HQET self energy was obtained in 
the three-loop one can be calculated using the methods of UTO). These calculations are 



based on integration by parts |11 



Using the identity (|^, we can obtain the result for the diagram of Fig. ||a without 
calculations: 



1 



Ca \ S(cj') - S(w) 
2Cl 



UJ' 



UJ 



(10) 



This result is also confirmed by direct calculation. Feynman integrals of the type of 
Fig. 0a, 

(11) 



{I ■ V + ujY^{i ■ V + uj'y^{pys 

can always be calculated, for integer z/i, z/2, by applying 



1 



{I ■ V + Uj){l ■ V + Uj') Uj' 



UJ 



1 



I ■ V + UJ I ■ V + uj' 



required number of times. We note that in integrals of the type (|rT|) with different 
velocities v have been examined. 
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k l + k k' 
Figure 8: Feynman integral V 



Calculation of Fig. requires more complicated Feynman integrals. A method of 
their calculation is presented in Sect. 0. Results for the HQET vertex, as well as various 
limiting discussed in Sect. ^. Sect. ^ contains a brief summary of the results 

obtained. In Appendix 0, we present general results for the Feynman integrals of the 
type of Fig. for arbitrary d and powers of the denominators. In Appendix ^ some 
issues related to the m — > oo limit of scalar integrals occurring in QCD are examined. In 
Appendix y, we discuss the relation between the QCD vertex at k,q <^ m and the HQET 
vertex. In Appendix 0, we present the HQET vertex for the heavy-quark scattering in 
an external gluon field in the background-field formalism. 



2 Triangle integrals 

2.1 Recurrence relations 

We consider the class of Feynman integrals (Fig. 



The cases i/q = 0, 1^2 = 0, i^i = are trivial, the results are proportional to 



V{0,1,1) = g{q' 



> ^2)d/2-2 r(2-c^/2)r^(c^/2-i) 



r(rf- 2) 



y(l,l,0)=XH, r(l,0,l)=J(a;')- 



(12) 



(13) 



When all the indices are non-zero, we use integration by parts |Tl|, similarly to [^0. 
Applying the operators {d/dl) ■ v, (d/dl) ■ I, and (d/dl) ■ {I — q) to the integrand of ([T^), 
we obtain the recurrence relations 



z/oO+ + 2z/il+(0" -uj) + 2z/22+(0" 





V 


= 0, 




V 


= 0, 


2-)' 


V 


= 0, 



(14) 
(15) 
(16) 
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where O^V^(z/o, z/i, z/2) = V^uq ± 1, z/i, 1/2), etc. When constructing the recurrence proce- 
dure we assume that all indices Ui are integer. 

If z/i < and z/2 7^ 1, we can raise z/i by (|1^); if z/i < and z/2 = 1, z/q 7^ 1, we can 
raise z/i or z/2 by (|I4D; if z/i < and z/2 = z/q = 1, we can raise z/i or z/q by (|16]). The 
case z/2 < is symmetric. If z/i > 1, we can lower it or z/2 by (p!6|); the case z/2 > 1 is 
symmetric. We are left with \^(z/o, 1, 1). 

Let us take times (p^if ), add times (|16D and tu' times (|15|), and subtract the 0~ 
shifted sum of (^5|) and ([T6|). We obtain at z/i = z/2 = 1 



Z/0O+ + (rf- 2z/o - 2) (cu + tu') - 2 (d- Z/Q - 2)0" 



^(^^0,1,1: 



l+2"(w-0") + 2+1" 







^(^^0,1,1) 



(17) 



where 



n = q"^ + Aiuuj' = (uj' + iuf - . (18) 

The integrals on the right-hand side of (|l^) are trivial. This relation allows us to raise 
or lower z/q. 

Therefore, all integrals (|T2D can be expressed, exactly at any d, as linear combina- 
tions of three trivial integrals (|T^) and a non-trivial one, 1^(1, 1, 1) = V{uj,uj',Q). An 
implementation of this algorithm in REDUCE can be obtained at 
http: / /wwwthep. physik.uni-mainz.de/Publications/progdata/mzth0101/ . 



2.2 Master integral 

The master integral 



V{uj,uj',Q) 



•^/^ J {l-v + u)P{l-qy 



(19) 



is convergent at c? = 4, except for the case = when it has a collinear divergence 
(infrared if g = 0). We shall consider the region < 0, cj' < 0, < 0, where no real 
intermediate states exist. 

Using the HQET version of Feynman parametrization (see, e.g., 0) 



a'^W T{a)T{l3) I {a + hyY+l^ 



1 T{a + (3) f y^-^dy 



(20) 



twice, we have 



where y and y' have the dimensionahty of inverse energy. 

For e = (rf = 4), calculating the integral in y' and substituting y = 1/z, we obtain 



7 log Az 

' (z + Q-cj + cj')(2-Q-cu + cu') ^ ' 







where 2; has the dimensionality of energy. Separating the logarithm as 

, -g2 - 2uj'z , z-2uo' , -q^ - 2uj'z 
.^.^o,.,^ = - log n / + log ■ 



z{z — 2u;) Q — — u;' z(Q — uj — uj') 

and making the substitution z = {—q^) / z' in the second integral, we obtain the repre- 
sentation 

J logjf^dz 

1^ 2 ' (J—U)—U]' 



^ {z + Q - UJ + uj'){z - Q - UJ + uj') 

7 log-^^d^ 
y (2 + g + cj-c<;')(-2-Q + cj-c<;')' 

which is explicitly symmetric under uj ^ uj'. The poles of the denominators at 2; = 
Q ± {uj — uj') are compensated by the corresponding zeros of the numerators. 
Finally, we obtain 

QV{uj,u',Q) =U2[- — +L12' * 



y 2uj J \ 2uj' ) 

, Q — UJ — uj' Q — uj + uj' 

+ log log 

^ -2uj ^ -2uj 

Q - UJ - uj' ^ Q + UJ - uj' tt^ ,^ ^, 

+ bg log . (24) 

^ -2uj' ^ -2uj' 3 ^ ' 

This expression has cuts at > 0, cj' > 0, and Q < \uj' — uj\., where real intermediate 
states exist. 



3 HQET quark-gluon vertex 
3.1 General results 

Here we present the one-loop HQET vertex, for arbitrary d and ^. The contribution of 
the diagram of Fig. ^ was given in (|10]). For Fig. ^d, we obtain 
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X \ 2{uj' + uj) [q^ + uJUj'iA + {d- 4)0j V(w, Q) 
q\A-i) + Auu;'{A+{d-A)i)\ g{q') 
+ {d-3)^{uj'I{uj) +ujI{uj')) 
9l i 



(25) 



|-2(rf-4)(cj' + w) 



X 



+ 



{d - Q)iQ'^q^uJuj' + nQ\2q^ + {A + {d- A)^)ujuj') 
+ {2 + {d- 3)0f^gWl V{uj, uj\ Q) 



A{d -3){d- 6)^gW(cj' + oof 

- 2Qq\q^ - Auju'id - A + {d - 3)0) 
+ nQ\{d-A){A+{d-A)0^ 

+ id~3)iA~id-A)0q')]Qiq') 
-{d-3){d- Q)^q\uj' + ujf{ujl{uj') + u'I{uj)) 
- {d-3)Q[2q\uj' ~uj){I{u}') -I{uj)) 

+ {A+{d- 4)0 (u;" - uj^){ujI{uj') - uj'I{uj)) 

+ 2{3 + Oq\ujI{u;') + u'I{u)) 



(26) 



where Q is defined in (|T8]). We have also derived the first contraction using (^ (Fig. §). 
The second contraction vanishes in the Feynman gauge, because the three-gluon vertex 
yields when contracted with v in all 3 indices. We also note that there are some 
cancellations in the generalized Yennie gauge, ^ = —2/{d — 3), as well as in the "singular" 
gauge ^ = —A/{d — A) (which was discussed in in connection with the three-gluon 
vertex) . 

When the bare vertex F'^ is expressed via the renormalized quantities 



9'o 



a 



-fx^'e^' [1 + 0(as)] , ao = a [1 + 0(«s 



(47r)^/2 47r 



it should become ZrT'^, where Z-p = 1 + Zia^/ {Atis) + ■ ■ ■ is a minimal renormalization 
constant, and the renormalized vertex is finite in the limit e ^ 0. Retaining only the 
pole parts Tioj) 2u/e, Q{q^) — ^ V{uj,uj',Q) 0, we obtain, either from ( P^ or 
from dl). 



[a-3)CF + ^CA 



as 
Aire 



(27) 
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.1/2 



When goT^ = gV^ Z^J'^ Zy is multiplied by the external leg renormalization factors ZqZj^ , 
it should give a finite matrix element. Using 



Za 
Zq 



1 - 

1 - (a - 3)Cf 



f) 



a - — ) + 3^^^^' 



[Zq follows from 



ni is the number of light fiavours) , we arrive at 

a. 



r7-2 ry -\ 



Po = y^A - -Tpni. 



(28) 



This means that the heavy-quark coupling with the gluon field in HQET is renormalized 
in the same way as the other QCD couplings. Of course, this must be the case, because 
otherwise renormalization would destroy gauge invariance of HQET. 



3.2 q parallel or orthogonal to v 

{uj' — uj)v, Q = 0. The denominators of 



In the parallel case q 
dependent. Inserting 



are linearly 



{l-qf + 2{uj' -uj){l-v + uj) 



into the integrand, we obtain 



1 



Uj' + UJ 



UJ — LO 



(29) 



exactly at any d. 

The vertex is, of course, proportional to f Therefore, we obtain, either from (p5|) 
or from (^), 



r. = 1 



Ca \ S(cu') - T.iuj) 



2C. 



uo' — UJ 



Ca 



9l 



(4vr)'i/2 A{uj' + uj) 



{d-d.) 



uo'liuo') — uoXiuo) 



(J - CJ)2 



UJ' — UJ 



Giiiu'-ujf 



(30) 



It has an imaginary part, which is contained in Q{{uj' — uj)"^). 

The case when q is orthogonal to v {uj' = uj, q^ = —Q^) does not lead to great sim- 
plifications. The contribution (|TD|) of Fig. |^a now contains dT,{uj) /duj = {d — 3)S(ti;)/a;. 
The contraction (^) is zero, and hence is parallel to v'^, too. The contribution of 
Fig. ^ to r„ is obtained from (^) by putting uj' = uj. 
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3.3 q on the light cone 

When = 0, the reduction algorithm of Sect. |2.1j breaks down. All V{vo,vi,V2) with 
z/q < vanish. We can use to lower z/q down to 1. Let us suppose that z/2 > 0; 
otherwise, we can interchange vi ^ 1/2, uj ^ uj' . The relation 



[{d — vq — ui — U2){uj' — oj) — uioj' + Z/2C1; + u;i^il^2 — uj'u22'^l ]V = 0, 



(31) 



which is uj' times (|T5[) minus uj times (]T6|), allows us to lower 1/2 down to 1. We are 
left with \^(l,z/i,l). Taking ([T6|), subtracting u' times (|l^) and adding 2uj' times l"*" 
shifted (0), we obtain at z/q = = 1 



- z/i - 3 + 2iu'{iu' -iu){d- 2z/i - 4)1+ 1^(1, z/i, 1) 
z/i - + 1)1^ 1^2-V{1, z/i, 1), 



(32) 



where the integrals on the right-hand side are trivial. Using this relation, we can lower 
or raise ui to 0. Therefore, all the integrals V{pq,ui,i'2) at = can be reduced to 
T{ijj) and Tioj'). For example, for V{1, 1, 1) we have 



(33) 



at arbitrary d. 

Repeating, with the new algorithm, the calculation of the diagram in Fig. ^ at 
= 0, we obtain 



9l 



1 



(471)^^/2 8(c/ - A)ujoj' 



X 



2(rf - 5 + (rf - 3)(d - A)^)uouj'{I{uj) - I{uo')) 
-{d-3){2 + {d- 4)0(c^"X(c^) - Lj^I{uj')) 
9l 



(34) 



Ca 



(47r)'^/2 I6{d - 6)wV2(tu' - uj) 



2{d - 6)uj^uj'\l{uj) -X{uj')) 
+ {2-{d- l)C)UJuj\uj''^I{uj) - uj'^I{uj')) 
+ {2 + {d- h)i){uj'^l{uj) - uj^X{uj')) . 



(35) 



These results can be also obtained from (|25|) , 
up to the g2 term. 



if we expand the numerator of (^) 
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The case g = belongs to all the categories considered above. We obtain, from each 



of the above results, = 0, 



3l (d - 3)^ [4 (2 + (d - 3)0 CF-{A+{d- 5)0 Ca] 



(36) 



{d - 3)J(u;) 



r„ = i- 

V(cj,cu,0) 
exactly at any d. 

3.4 VL = {) 

Another interesting case is f2 = (g^ = —Aujoj'). After reducing V^(t'o, J^i, ^^2) to V^(z/o, 1, 1) 
and trivial integrals (Sect. we can use (|1^ to reduce to 0. In particular. 



V{ijJ, ijj' , ijj' + a;) 



rf- 3 



I{uo) I{uo') 



2{d - A){to' + uj) [ ~ ■ UJ 

for any d. Repeating the calculation of the vertex, we obtain 

1 



UJ' 



9l 



X 4 



+ 



+ 



(47r)rf/2 i6(d - 4)(d - Q)uuj'{uj' + 
-{d-Q){A + {d-A)i){uj' + ujf 
+ 4(c/ - 3)(d - 4)^cjcj']cucu'(cj' - cj)^(-4cjcu') 
"(4(c/ - 6) + (d - 3)(3rf - 16)0(^' + 

- A{d - ?,){d - A)^ujuj' 

x{d- A){uj' + uj){uj'I{lo) - ujI{uj')) 
{A{d - 6) - (rf - 3)(c/ - A){d - 8)0iu}' + 

- A{d - 3){d - A)^ujuj' 



X {uj' - uj){uj'I{uj) + }, 



b A ^(47r)'^/2 32(d-6)(rf-8)cjV2(cu' + cu)2 



{d -6){d- 8)(4 + (d - 4)0(cu' + uj)^ 
+ A{d - 8){2d -{d-3){d- A)^)ujuj'{uj' + ujf 



(37) 



(38) 
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+ 



+ 32{d - 3)(rf - 8 - 30uj^u'^]uju'g{-4uuj') 
{2{d - 7){d - 8) - {d - A){d - 5)0(^' + ujf 
-A{d-5){d-8-3^)ujuj' 

x{d- 3){uj' + uj){uj'^I{uj) + uj'^I{uj')) 
{d - 5){2{d - 8) + (rf - 4)0(cu' + 

+ 4(d - 8 - 3^)uju' 



(39) 



3.5 Quark(s) on the mass shell 

If one of the quarks is on its mass shell, say uj' = 0, then 

gv(.;,o,g) = 2Li,(;^j+iog ^-^log qT^- — 



(40) 



at (i = 4. The contractions of the vertex are obtained from (|2|), (|2|) by setting J(c<;') = 
and then u' = 0. 

When both quarks are on shell {u = u' = 0), 



r^g^ = 0, 



i[-2(4rf-13) + (rf-4)e]^(g^). 



(41) 



In this case, V(0, 0, Q) does not appear. 

It is easy to consider the cases when q is parallel to v and u' = 0, and when = 0, 
uj' = 0. 



4 Conclusion 

We have obtained general expressions (p!OD, (p5|) and (p6D for the one-loop HQET quark- 
gluon vertex. Using recurrence relations ([T^ ) -([T7| ) , we expressed the results in terms 
of one non-trivial integral V{uj,u',Q) (p!9| ) and some trivial integrals (|13]). For the in- 
tegral (|T9|) in four dimensions, we have obtained an analytic result (|2^) in terms of 
dilogarithms. In Sections |3]^-|3]^ we have also studied some special limits of interest. 

In Appendix ^ we have provided some results for the integrals (p!2[) with arbitrary 
indices and in arbitrary dimension. We have also discussed, in Appendix ^, how the 
HQET limit can be obtained directly from the standard integrals occurring in the QCD 
calculation. Using this prescription, in Appendix |C| we have examined the m — > oo limit 
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of the general QCD result for the one-loop quark-gluon function, and we have found 
that it is in agreement with our calculation. We have also presented the result for the 
background field vertex (Appendix 0). 
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A One-loop HQET integrals 

Result for the two-point HQET integral is well known 0, 

i r d'^l 



7 {i-v + ujy'->{py 



, .y^+u^u,, o xd-.o^2. r(z/o + 2z/-rf)r(rf/2-Z/) 



For the triangle integral (|12[), using Feynman parameters, we arrive at the following 
double integral representation: 

r(z/o) r(z/i) r(z/2) 



oo oo 



^ dy dy' y''--' y'"'-' jy + y'yo+'^^+-^-<i 







[1 - 2ujy - 2u'y' - q^yy 



The symmetry (cu, z/i) ^ {u', U2) is explicit. 

In the special case a; = a;' = 0, the integral (|^) can be evaluated in terms of F 
functions, 

r(^o/2)r(^/2 - z/0/2 - z/i)r(rf/2 - ^0/2 - v^) 

r(z/o)r(z/i)r(z/2)r(rf-z/o-z/i-z/2) ' ^ ^ 

In particular, for z/j = 1 we get 

V(O.O.Q) ^ - ^'^^^'"»,^-_^'3^,^:r- ^>''^> ^ + (45) 

In another special case, = 0, we get a hypergeometric function, 
V^(z/o,z/i, z/2) 1^2=0 = (-l)-o+.i+.22-o(_2a;)'^-'^«-2'^^-2^^ 
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r(z/o + 2z/i + 2z/2 - d) r{n/2 - ui - U2) 

X 



X 2F1 



r(z/o) r(z/i + U2) 



Vl + Z^2 



UJ 



1-- . (46) 



Using simple transformations of 2F1 function, it is easy to see that the result obeys the 
symmetry {uj.i'i) ^ (cu', 1/2), as it should. Note that the structure of the result (|i6| ) is 
quite similar to that of the two-point integral with different masses and zero external 
momentum, see eq. (2.9) of ||l5l. When 1/1 = 1/2 = 1, the result (^61) reduces to 



V{uo, 1, 1)1,2=0 = {-ir2^~'r{uo -d + 3) T{d/2 - 2) 



(47) 



Let us represent the denominator of ( ^31) in terms of double Mellin-Barnes integral, 
expanding with respect to u and u' (see, e.g., Eq. (3.4) of |jl5|). Then, the resulting 
momentum integral can be recognized as 

V{uo + tl + t2, 1^1 + h, V2 + t2)L=^,=o, d^d+2ti+2t2' W 

where t\ and ^2 are the contour integration variables. Using ( ^4]) we can evaluate the 
integral (^Sp in terms of F functions. Making a linear substitution for the contour 
integration variables (ti = s + 1, ^2 = s — t), we arrive at the following double Mellin- 
Barnes representation for the integral (0): 



r(z/o) r(z/i) r(z/2) v{d -v^-v^- 1^2) 



100 100 / i\ s 

is it -ii^ (^) ri-s-t)r(t-s) 



(27ri)2 J J \ g2 

—ioo — ioo 

xT{d/2 - uo/2 -ui-t) T{d/2 - uo/2 - z/2 + t) 

xr(z/o/2 + z/i + z/2 - d/2 + s) r(z/o/2 + s) . (49) 



B QCD integrals in the HQET limit 

Here we discuss the relation of massive integrals occurring in standard QCD calculations 
and HQET integrals. First, let us consider the two-point integral with one massive line 
and one massless line. 



For general values of z/q, z/ and d, such integrals have been examined in ||16|,[L 
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When we substitute p = mv + k, the massive denominator in ( |50| ) becomes 



2m(/ + k)-v + {l + kf 



(51) 



For k <^m, several regions of integration in / are essential. When / ~ fc, we can expand 
the heavy propagator (^) in both k/m and l/m, and it becomes the HQET propagator. 
The leading term of this HQET contribution (called "ultrasoft" in |]l8| , p!9[] ) yields l/rrf'^ 
times an HQET integral (^21) , which is proportional to {—2ujY~^°~'^'^ , by dimensionality. 
Higher terms form an expansion in k/m. 

Let us subtract and add this expansion of the heavy propagator to the exact one. In 
the difference, the contribution of small / ~ A; is suppressed; typically, / ~ m. Therefore, 
we can expand this integrand difference in regular series in k/m, and integrate term by 
term. Integrals of all terms of the HQET integrand expansion in k/m vanish in dimen- 
sional regularization, because they contain no scale. Therefore, this "hard" contribution 
can be obtained by expanding the exact QCD integrand (^T]) in k/m, and integrating 
term by term. It is analytical at /c = 0, by construction. The leading term is proportional 
to m'^~'^^°~'^'^ , by dimensionality, whereas the higher terms form an expansion in k/m. 



This separation of J{h'Q, v] m, 0) at A; -C m into two contributions |]T^ is a particular case 
of a more general threshold expansion 



Note that k plays a role of the threshold 
parameter, since — m^ ~ mk. 

We can check these qualitative considerations, using an explicit expression for J(z/o, z/; m. 
It can be presented in terms of 2-F1 function oi z = /m^ (see Eq. (10) of ||16|)- Note 
that in the HQET limit z approaches 1, 



[mv + fc)" 



2uj fc2 



m^ 



m m 



2 ' 



i.e., it is at the border of convergence of the 2F1 function. Transforming from the variable 
2; to 1 — 2;, we obtain (see also Eqs. (1.12)-(1.15) of [l] 



J(z/o,z/;m,0) = ivr^/2(_i)-o+-^<i-2vo-2^ 

f r(z/o + v- d/2)T{d -vq- 2u) 



[ r(z/o)r(d - z/Q - z/) 

Y(uo + 2v- d)T{d/2 - v) 



2F1 



z/, z/q + z^ — d/2 



vq + 2v -d + l 



\-Z 



X2F1 



d/2 — V, d — Uq — u 
d — vq — 2v + 1 



1 - z 



(52) 



We see that the first term here is nothing but the "hard" contribution. It has a 
prefactor m'^"^'^""^''. The prefactor of the second 2-^1 function is 



m 



d—2vo—2v. 



1~Z) 



m 



-2tu) 



d~vo—2v 
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up to higher powers of 1/m. This is the HQET ( "ultrasoft" ) contribution; in the leading 
order, it yields 

J(z/o, z/; m, 0) m'^''^{2m)-''''I{pQ, z/), 

where liyo^p) is defined in (^). 

The value of {2mY° J{vQ,v]m,Q) at the singular point 1/m = is (up to a factor 
ivr''/^) the HQET integral When setting 1/m = 0, we discard the "hard" contribu- 
tion, which is proportional to m'^"^""^'", because we can always choose d small enough for 
this contribution to vanish. The naive Taylor expansion in 1/m is given by the HQET 
contribution. Similarly, the value of J(z/o, z^; 0) at the singular point = is the 
QCD on-shell integral. When setting = 0, we discard the HQET contribution, which 
is proportional to (— 2a;)'^~''°~^^, because we can always choose d large enough for this 
contribution to vanish. The naive Taylor expansion in k is given by the "hard" contri- 
bution. Neither of these naive expansions, taken separately, describes J(z/o, z/; m, 0) at 
k <^ I. It is given by the sum of both contributions, "hard" and HQET ("ultrasoft") 
ones (see in |T^). 

Now, let us consider off-shell three-point integrals with one massive and two massless 
internal lines. For arbitrary powers of propagators, such vertex integrals has been con- 
sidered in |T|] (see Eqs. (25)-(29)) and (Eqs. (4.1)-(4.5)). Using the Mellin-Barnes 
representation, Eq. (4.2) of (see also Eq. (26) of |[l^), we see that two (of three) 
arguments, pl/m^ and p^/m^ (we denote p\ = ^23, P2 = ^13? P3 = ^12) ^I's approaching 
1 in our limit, which is at the border of convergence of the corresponding function. 

Again, as in the two-point case, we need to construct analytic continuation of this 
representation, in terms of the variables containing (m^ — pi) and (m^ — P2). To do 
this inside the Mellin-Barnes representation, we use the contour-integral version of the 
corresponding formula for 2-^1, namely, Eq. (A7) of (where we put z/3 = z/q). In this 
way, we arrive at 



r(z/i)r(z/2)r(z/o)r(c?- z/i - z/2 - z/q) 



1 [Tf, , , frn^ — p^y/rn^ — plVf ply 



{27ri)^JJJ \ m? J\ m? J\ m? , 

— ioo 

xr(-s)r(-t)r(-M)r(rf-z/o-2z/i-2z/2-s-t-2M) 

xr(z/o + V1 + U2- d/2 + s + t + u) 

xV{ui+t + u)V{v2 + s + u) . (53) 

Now, we need to analyse the contributions of the poles in the right half-plane of the 
contour variable u, since they correspond to increasing powers of p\jrr? . There are only 
two series of poles (with j = 0, 1, 2, . . .), 

(i) u = j (due to r(— n)) and 
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(ii) u = dl2- z/o/2 - ui - U2 - s/2 - t/2 + j/2 

(due to V{d — uq — 2ui — 2u2 — s — t — 2u)). 

The prefactors contain TTi'^-^j^o -21^1 -2^2 (^c^gg (^j)) and m"'''^ (case (ii)). The series (i) yields 
the "hard" contribution. The HQET ( "ultrasoft" ) contribution is given by the series (ii), 
whose leading term is j = 0. These two contributions are related to the singular limits 
1/m = and A; = g = in a way similar to the two-point case. 

Picking up this leading HQET contribution, we arrive at a double Mellin-Barnes 
representation (in terms of the remaining contour integrals over s and t), where the 
kinematical variables involved yield (in the limit m 00) the HQET variables, 

— p\ —2uj' — P2 —2uj 

After introducing the new contour variables (s + t)/2 and (s — t)/2 we see that the 
resulting Mellin-Barnes representation is equivalent to (^91), so that 

Ji(z/i, z/2, z/q; m) ^ i7r"'/^(2m)"''°V(z/o, z/i, z/2) . 

Therefore, the HQET integrals can be formally obtained directly from the standard 
loop integrals, by picking up the formal Taylor series in 1/m which has no prefactors 
containing m to the power depending on d. Such prescription is similar to some other 
prescriptions in dimensional regularization. For instance, considering the massless limit 
of the integral (|50|) we need to represent the result in terms of the functions of the variable 
w? /p^, and then discard the contribution containing m~'^^ (see Eq. (11) of [|16|). We have 
also demonstrated that the HQET contributions are equivalent to the "ultrasoft" ones, in 
the language of the threshold expansion [|18|. In other words, in the cases considered the 
exact result is given by the sum of two contributions. The first one is given by a formal 
Taylor expansion of the integrand in the small parameter of the threshold expansion (the 
"hard" contribution). The second one is nothing but the HQET series in 1/m. 

C QCD vertex at m ^ 00 

If we substitute the "hard" parts of all scalar integrals into the QCD vertex, then, in the 
limit A; — > 0, g — ^ 0, we just get the on-shell vertex at g = 0. Corrections to this limit are 
regular expansion terms in /c/m, q/m. Since we do not consider 1/m suppressed terms 
here, these corrections can be omitted. 

If we substitute the HQET parts of all scalar integrals (see Appendix |BD, we should 
obtain the HQET vertex, which was calculated in Sect. ^. In order to make a strong 
check of both the results of 0] (where the one-loop quark-gluon vertex was calculated 
in arbitrary gauge and dimension) and of the present ones, we consider here the HQET 
limit of the QCD vertex [|]. 
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Using the standard decomposition of the quark-gluon vertex (see also in [P,p2||;, 
it can be spht into longitudinal and transverse parts, 



r'^ = EA.^r + E^^^'> (54) 

i=l 1=1 

where Aj and r, are scalar functions depending on kinematical variables, whereas and 
T/" are vectors which also involve Dirac matrices (see Sect. IID of for further details). 
We substitute 

Pi = —mv — k — q, P2 = fnv + k, Ps = q. 

At the leading order in 1/m, the initial and final quark spinors obey ■^u = u. Therefore, 
we can sandwich between the projectors {'^ + l)/2, and use 

2^ 2 2 2 ■ 

With the required accuracy, there are three independent structures: 

L'^^v^' + 0{l/m), 

n - {uj' - uj)q^ + 0(l/m), Ti: ^[Y, i]- 

The others are 

^4m2L^ + 0(m), ^ -2mL^ + 0(1), ^ C(l), 
Tf ^ mT^ + Oil), 2m^T^ + C(m), 

T^^O{m), T^^O{l), Tj^^O{m), ^ niT^ + 0{1). 

According to Appendix 0, the HQET-relevant "ultrasoft" parts of the scalar integrals 
are 



liuo) I{uo') 



2\ 



2m ' 2m 



Wi 7, > W2 



2m ' (2m) 2 uj' — uj ' 

where the notations y^j, Ki^i and rj are defined in Sect. IIA of [^. 

Therefore, at the leading order (l/m°) the HQET limit of the QCD vertex (|5^) yields 

(Ai + 4m2A2 - 2m\^)Lt + (mn + 2mV2 + T3)T.^ + (rs + mT^)T^ . 



19 



Using the results for Aj and x; listed in , we have obtained that at the order 1 / mP the 
coefficient of the chromomagnetic structure TJ^ vanishes, whereas those of Li and 
reproduce the results ([T0|), (|25|) and (^), for arbitrary d and ^. 

Thus the QCD vertex at small k, q is equal to its on-shell value at = g = plus 
the HQET vertex, up to 1/m corrections. We can reformulate this statement: the QCD 
vertex in the on-shell renormalization scheme is equal to the HQET vertex in the on-shell 
renormalization scheme, up to 1/m corrections. In QCD, the on-shell renormalization 
subtracts from the one-loop correction its value at /c = g = 0. In HQET, the on-shell 
renormalized vertex equals the bare one, because its value at k = q = vanishes. 



D Background-field vertex 

The background field formalism is convenient for considering heavy quark scattering 
in an external gluon field. In this method, the vertex (Fig. |lD with the background gluon 
differs from the ordinary one. The diagram of Fig. ^ still gives (p!0D , because the quark- 
gluon elementary vertex does not change. The contraction of the three-gluon vertex 
with a background gluon momentum contains no ghost terms, just the first difference in 
Fig. Therefore, 



Ca 
2Cf 



This is also confirmed by a direct calculation. For the other contraction, we obtain 

1 



X \ 2{uj' + Lu) [2fi2g2 _ _ 8)^^]2g2 _ ^^^q2^2 

- (16 + {d^ - 18d + 40)0 
-{d-3){d- A)('nq^ujuj' - d^i'^Q^q^uJuj' 
-8{5d- 12)^2g2^ V^] V{uj, uj', Q) 

m^q^ + (4(2rf - 7) - (rf - 4)0e^^Q^ 

- 8(rf - 3)(2 - (rf - 5)0^fi<5W 

+ 4(d-3)(rf-4)eWw^' 

- 16(rf - 3)(rf - 6)e'QVa;'']e?(g2) 

2(rf - 3)^ [(2 + i)q^{uj'^I{u) + u^I{u')) 
+ 2(4 + (rf - A)i)uju}\uj'^ - uj'^){uj'X{ij) - 

- ((1 + 0^ + {d- Q)^ujuj')q\uj'l{uj) + ujI{uj')) 
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- {<f + (6 + {d-^)i)uju')q^{uX{u) + uj'I{u'))\ | 
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